Estimating the severity of landscape degradation in future management scenarios based on modeling the dynamics of Hoor Al-Azim International Wetland in Iran-Iraq border

Temporal and spatial changes in land cover in wetland ecosystems indicate the severity of degradation. Understanding such processes in the past, present, and future might be necessary for managing any type of development plan. Therefore, this research has monitored and analyzed the Hoor Al-Azim International Wetland to determine the orientation of its changes in various future scenarios. Wetland status modeling was conducted using developed hybrid approaches and cellular automata along with evaluating the accuracy of the modeled maps. The dynamics of the landscape were simulated using a higher accuracy approach in three scenarios—Water Conservation, Water Decreasing, and Business-as-Usual- to get the level of degradation of the wetland. The results showed that the amount of water in the wetland has decreased in all three periods, and the salt lands and vegetation have undergone drastic changes. The water bodies experienced a reduction of 148,139 ha between 1985 and 2000, followed by a decrease of 9107 ha during the 2000–2015 period. However, based on the results, these developments are expressed better by the developed hybrid approach than the CA-MC approach and are more reliable for future simulation. The figure of merit index, which assesses the hybrid model's accuracy, yielded a value of 18.12%, while the CA-MC model's accuracy was estimated at 14.42%. The assessment of degradation in hexagonal units showed the least degradation in the water conservation scenario compared with the other two scenarios in 2030.


Study area
Hoor Al-Azim wetland is a part of one of the largest wetlands in West Asia, called the Mesopotamian wetland, with an area of over 9000 km 2 .This wetland (Hor, marsh) is composed of the central Hor, Horalhammar (Hammar Marshes), and Horalhawizeh (Hawizeh Marshes).Central Hor and Horalhammar are entirely located in Iraq, and about two-thirds of Hawizeh is in Iraq and one-third is in Iran.The Iranian part is known as Hoor Al-Azim (Hadith 35,36 ).Hoor Al-Azim wetland, with an area of 771,161 ha, is in the southwest of Khuzestan province, in the northwest of Bostan city, and near Chazabe village at the coordinates of 71° 51′ E and 71° 76′ W (Fig. 1).The Hoor Al-Azim wetland is a part of the permanent freshwater wetland in lower Mesopotamia, which is located in the borderline region between Iran and Iraq.The wetlands of Mesopotamia are one of the 200 most important biological areas in the world 5 .The Hoor Al-Azim wetland is important internationally for historical, cultural, and environmental reasons.However, it has always been directly and indirectly affected by anthropogenic degradation at the national and international levels.The primary threats to the water ecosystem from human activities, including dam construction in Turkey, extensive well drilling in Iraq, and the Karkheh Dam in Iran, all of which have significantly reduced the water levels of this vast ecosystem.Additionally, the degradation of the Hoor Al-Azim wetland is exacerbated by oil exploration, road and dyke construction, river deviation, and the aftermath of wars (Iran-Iraq and Gulf wars), alongside natural factors like recurring droughts 37 .The reduction in wetland area has directly contributed to vegetation loss and surface diminishment, leading to increased susceptibility to dust particles and the generation of fine dust 38 .

Data preparation and mapping of modeling variables
Landsat satellite images from 1985, 2000, and 2015 produced by Nadaf 39 were used to simulate land cover change dynamics in the Hoor Al-Azim wetland.During the pre-processing phase of satellite images, geometric, atmospheric, and radiometric error corrections were conducted.The fuzzy classification was utilized in the eCognition software to generate maps.Each classification class in eCognition is planned based on the class descriptions, and each description includes a set of fuzzy statements that allow the assessment of areas and logical functions.The output of this system has two parts as follows: 1.The fuzzy classification with partial information of classes provides a fragile classification for placing objects in a class.2. The fuzzy principle can have individual conditions or a combination of several layers that include all evaluations for classes.In eCognition, conditions are defined by statements that are entered for class descriptions.Definitions related to individual conditions, logical combinations of multiple conditions, and nearest neighbors are given below 39 .
(a) Individual conditions: Individual conditions are defined by a one-dimensional membership function.For visual object features, the one-dimensional membership function is defined by the graphical surface.Through the one-dimensional membership function, all available information about the relationship between areas and the assignment of classes can be completed.
(b) Combined conditions: in most cases, the principles required for land use classes cannot be described through individual conditions.Classes are often composed of combined conditions that are determined by the operator.In the area beyond, individual conditions can be expanded.
(c) Nearest Neighbor: Fuzzy perception from the point of view of the nearest neighbor, which is used in eCognition, automatically generates the multidimensional membership function.Five land use classes in the study area include humid land, water body, salt land, sparse vegetation, and dense vegetation (Table 1) (Nadaf  et al. 2015).
After assigning suitable weights to the image bands, segmentation was conducted.During this process, a scale parameter of 10 was selected, with homogeneity preferences set at 0.8 for color and 0.2 for shape, alongside 0.9 for shape smoothness and 0.1 for compactness.
The accuracy of the classified land cover maps was assessed using indices derived from the error matrix 40 .To enhance the comprehensiveness of the study area's coverage, ground truth points along with their land cover data were meticulously gathered.Each sample site was precisely marked with GPS coordinates and corroborated with ground truth information.The accuracy of the classified maps was rigorously evaluated by comparing the results of classified map data against the ground truth on a per-category basis, utilizing an error matrix to delineate the classification results' reliability.Accuracy assessment of the classified maps was done based on the error matrix and Kappa Statistic with about 250 random samples of ground truth data for each year 41 .Kappa was calculated using Eqs.(1) 42 .
where r is no. of rows in the error matrix, X ii is no. of observations in row i and column i, X i+ is total observation in row i and X +i is total observation in column i.

Variables
Slope, digital elevation model, distance from main roads, distance from the river, distance from sparse vegetation, distance from dense vegetation, distance from humid lands, distance from water bodies, distance from salt lands, evidence likelihood to change the map, annual average temperature, and precipitation variables (Supplementary Table 1) were produced to select effective variables in the expression of changes in Hoor Al-Azim wetland (Supplementary Fig. 1).Topography has been considered by other studies as an effective driver variable determining land use and land cover changes 43,44 .Distance from water bodies and vegetation were important factors in modeling changes in the Saudi Arabian eastern coastal city of Dammam 45 .Distance to the rivers was recognized as an essential underlying driving force of land cover change in Iranian Northern Zagros forests 46 .Proximity to main roads and evidence likelihood to change map have been reported as essential variables in modeling the dynamics of land cover in the Zanjan Province, Iran 19 .The total annual precipitation is considered an explanatory variable in the generation of land use scenarios in the Catamayo Chira basin, located between Ecuador and Perú 11 .
The digital elevation model of the studied area, with a spatial resolution of 30 m, was extracted from the ASTER digital elevation model.The slope layer was also produced using a digital elevation model map.Sparse vegetation, dense vegetation, and rangeland were extracted from the land cover image of the older year (in the model input) as Boolean layers, and the variable representing the distance from each of these classes was prepared 47 .The evidence likelihood to change map was generated using the evidence likelihood transformation tool.In this way, the map of changes during the study period was produced as a Boolean layer and placed in the Transition or land cover layer section.Also, the land cover map of 1985 (a map of an older year in the model's input) was entered into the variable input section of the transformation panel, and the evidence likelihood to change map was generated.

Research method
This research was done in four steps, according to Fig. 2. www.nature.com/scientificreports/ Step 1: Monitoring and analysis of land cover changes The analysis of changes in the study area was calculated using LCM for the first period (1985-2000) and the second period (2000-2015).LCM can determine the change process based on location, calculate transition probabilities between land cover classes, and predict future land cover maps 44 .These changes include gains, losses, and net changes for each land cover class and the transition from one class to another.Also, maps of the spatial trend of changes from water bodies to other land cover in Hoor Al-Azim International Wetland were produced using the results obtained from the analysis of changes.
Step 2: Modeling of the wetland status

Hybrid approach
The present study integrates the approaches of the artificial neural network, Markov chain, and multi-objective optimization into a single model.LCM uses transition potential maps for each transition sub-model of land cover classes to assign simulated transitions based on land cover changes in the past 46 .In the present study, a multi-layer perceptron artificial neural network was used to model the transition potential.The network uses feed-forward algorithms to calculate weights for input values, input layer nodes, hidden layer nodes, and output layer nodes, which are transferred to the hidden layer and the output layer 48 .In this study, the number of repetitions is set at 10,000 to achieve high accuracy.Accuracy rate statistics and Skill statistics are also used to evaluate the neural network.The value of the Momentum Factor in the present study was considered to be 0.5.The LCM model uses Cramer's V to examine the relationship between each variable, the distribution of land cover classes, and the changes that have occurred.The value of this coefficient is between 0 (no correlation) and 1 (complete correlation) 49 , which is calculated using Eq.(2).
where X 2 is the chi-square statistic, n is the number of samples, and k and r are the number of rows and columns in the probability table, respectively 47 .After determining the power of the explanatory variables using the Cramer coefficient (a high Cramer coefficient means that the explanatory potential of the variable is good), ten variables were considered, and the artificial neural network with slope variables, digital elevation model, distance from main roads, distance from the river, distance from sparse vegetation, distance from dense vegetation, distance from humid lands, distance from water bodies, distance from salt lands, and the evidence likelihood to change map was instructed.
In order to ensure the selection of the best sub-models for modeling, a spatial analysis was performed on the change pixels using explanatory variables of the changes 50 .The sub-models used in the research include the water body to sparse vegetation, water to dense vegetation, water to humid land, water to salt land, humid land to sparse vegetation, dense vegetation to sparse vegetation, dense vegetation to humid lands, and salt land to humid lands.
The results of the artificial neural network are combined with Markov chain and multi-objective land allocation to simulate the spatio-temporal dynamics of land cover changes in the Hoor Al-Azim wetland.In this study, the Markov chain model was used to estimate the rate of change that can occur in a specific future year based on previous land cover maps and the probability of change from one land cover class to another.A prediction map for 2015 was created using the Markov chain to evaluate the accuracy of the model 51 .After ensuring the accuracy of the model, the land cover maps for the years 2000 and 2015 were used to predict the land cover maps for 2030 44 .

CA-MC approach
The CA-MC model is a combination of cellular automata, Markov chain, and multi-objective land allocation, which are used to predict future changes in land use and land cover 52,53 .In this research, the CA-MC approach was used in 1985 and 2000 to predict land cover changes in 2015.The basic assumption of the approach is that the land cover status in the future (t + 1) can be defined as a function of the current land cover (t) as shown in Eq. ( 3), where x (t+1) shows land cover at time t + 1 and x shows land cover at time t 54 .The structure of the Markov chain model for modeling land cover changes includes a vector (N) with dimensions Mx1 (M represents the number of land cover classes) that describes the current land cover distribution and a MxM matrix of transition probability (P) between each pair of land cover classes i and j.The transition probability is calculated according to Eq. ( 4).

Step3: Accuracy assessment and future scenario planning
The prediction error and accuracy of the model were evaluated based on three land cover maps: the actual maps from 2000, and 2015, and the map generated from the 2015 model simulations.The FOM index was calculated according to Eq. ( 5) to evaluate its accuracy.The value of this index is between 0 and 100, where 100 indicates a complete overlap between the ground truth map and the simulated map 55 . where

Scenario creation
Three scenarios were considered for modeling land cover changes in 2030, and the labels and selected sub-models for each are described below.
Scenario 1 (Water Conservation) According to the sixth development plan (2016-2021) 56 in terms of social, economic, and environmental sustainability in Iran, we selected the sub-models in order to increase the transition probability of other classes to water.In this scenario, the sub-models used include sparse vegetation to water, dense vegetation to water, humid lands to water, and salt lands to water.

Scenario 2 (Water Decreasing)
There is a risk of the wetland drying up due to upstream dams of the wetland, oil facilities construction, new roads development, and water exploitation 57 .This scenario assumes high economic growth and predicts continued shifts in consumption patterns.In this scenario, the sub-models used include water body to sparse vegetation, water to dense vegetation, water to humid lands, and water to salt lands.
Scenario 3 (Business-as-Usual) This scenario assumes current trends and projects the wetland's future status for 2030.It is selected to explore the sensitivity to recent drivers of changes using sub-models such as water to sparse vegetation, dense vegetation, humid lands, and salt lands; humid land to sparse vegetation; dense vegetation to sparse vegetation, humid land; and salt lands to humid lands.

Step 4: Evaluation of the landscape degradation
To assess the degradation of the landscape, the study area was divided into a network of uniform hexagonal units.Hexagons are the closest regular geometric shapes to the circle, which can cover an area with no overlap and have a smaller perimeter than a square with the same area, so it reduces the error caused by the edge effect 58 .In this study, due to the large scope of the study 59 and the quality and accuracy of the data 60 , hexagons with an area of 40,000 hectares were selected as management and monitoring units after numerous tests and visual checks to ensure that all changes were included.Management and monitoring units should be large enough to represent different types of land cover 61 and to reflect the severity of degradation caused by land cover change and landscape metrics in temporal and spatial gradients.The quantification of CA, MPS, ED, MPFD, and MSI metrics www.nature.com/scientificreports/ was conducted to evaluate the degradation of the landscape in the study area using the Patch Analyst tool, and the details of the metrics are presented in Table 2.The application of landscape metrics and the scientific basis of their relationship with land degradation is a precedent in various studies 32,30,31,59,62 .In this research, the severity of degradation is a relative value obtained from the value of the metrics that influence landscape degradation.These metrics are a measurement category of interval scales have converted into ordinal scales during mathematical processing.In each of the hexagon units, the rate of degradation was calculated based on Eq. ( 6).
where KI determines the intensity of activities in the study unit, K represents the degradation intensity class of each metric, I is the metric, and V i represents the intensity of ecological vulnerability.In this research, the value of K, which represents the degradation intensity class of each metric, has been obtained after calculating the landscape metrics and normalizing them based on the fuzzy logic in Eq. ( 7).
In Eq. ( 7), Z i represents the normalized value of the metric in the unit of hexagon i. X i is the metric quantity in hexagon units i.The min(x) indicates the minimum metric quantity between all the units of the hexagonal network.Degradation severity codes were assigned based on the impact of suitable and unsuitable land covers on the wetland landscape and the effect of metrics, as shown in Table 3.
The land cover degradation index was qualitatively determined between 0 and 1 by expert opinion.Lower values of this index indicate less degradation, and vice versa (water: 0, humid lands: 0.1, dense vegetation: 0.2, sparse vegetation: 0.7, and salt lands: 1).Vulnerability refers to the degree to which a system, subsystem, or system components are susceptible, unable to cope with its adverse effects, and are damaged when exposed to stimulating factors (disturbance or pressure) 64 .The intensity of ecological vulnerability was also calculated based on the percentage of suitable land cover in each hexagonal unit.

Analysis of changes in land covers
The dataset used for modeling the wetland landscape dynamics includes variables and land cover maps from the years 1985, 2000, and 2015.These maps were meticulously prepared using GIS technology and are illustrated in Supplementary Figs. 1 and 2.An overall Kappa coefficient of 0.8546, 0.9110, and 0.8887 was obtained for the Landsat images of 1985, 2000, and 2015, respectively.The detailed results of the confusion error matrix, including errors of commission, omission, and Kappa Statistics in each class, are presented in Supplementary Table 2.

Group metrics 2 (ED, MPFD, MSI)
The water bodies accounted for 25.1% of the total area in 1985.By 2000, this percentage has decreased by 7.1%, resulting in an annual reduction rate of 8.3%.The gain and loss in water bodies were 15,735 ha and − 163,487 ha, respectively (Fig. 3).From 2000 to 2015, the net change from dense vegetation, water, and salt lands to sparse vegetation was 18,034, 9985, and 62,946 ha, respectively.The area of water bodies shrank by 9,107 hectares, corresponding to an annual reduction rate of 1.1%.Sparse vegetation increased during this period at an annual rate of 2.5% and 91,509 ha (Supplementary Table 4).The net change from water to humid and salt lands was 9,741 and 7,449 ha, respectively.According to the spatial trend map (Supplementary Fig. 3), the most significant changes from water to salt lands occurred in the southern region during 1985-2000 and in the centeral and southern parts of the studied region during 2000-2015.

Modeling the transition potential and hybrid approach
The importance of the variable in the artificial neural network model was measured using Cramer's coefficient, which ranges from 0 to 1.According to Cramer's coefficients in Table 4, all variables showed "good" power in describing changes.According to Supplementary Table (3), the water to salt lands sub-model showed the highest accuracy rate (98.04%), and the dense to sparse vegetation sub-model showed the lowest accuracy rate (82.22%).According to the results of the Markov chain, the water with the highest probability will turn into humid lands from 2015 to 2030 (0.3633).The map modeled by the developed hybrid approach is depicted in Fig. 4a.

Modeling changes with the CA-MC approach
The changes in land cover area were obtained using Markov chain analysis between 1985 and 2000, based on the area transition matrix in Table 5.Finally, the CA-MC approach was implemented using the area transition matrix.Land cover maps were repeated 10 iterations with a 5 × 5 filter, and the land cover change map for 2015 was depicted as Fig. 4b.

Assessing the accuracy of modeled maps
The FOM index range, derived from evaluating the accuracy of the hybrid model, is 15.6%.This value signifies the high efficiency and capability of the model utilized in this wetland, attributed to its substantial size and the varied changes during the calibration period in contrast to the simulation period.Also, the accuracy of the CA-MC model was estimated to be 14.42% (Fig. 4).The Hits, Misses, and False Alarms zones are shown in Fig. 4.

Scenario planning of the future status of the wetland
After evaluating the accuracy of the model, different land cover scenarios for 2030 were produced based on the three policies outlined in the methodology, utilizing the superior hybrid approach identified (Fig. 5).

Evaluation of the landscape degradation
In Fig. 6, the hexagonal network consisting of 25 units is depicated as the basis for evaluating degradation in the wetland.This map was created by overlaying the land cover map with the hexagonal network.Additionally, a desired network was developed for future scenarios.The rate of landscape degradation was calculated using to the mentioned method based on landscape metrics in Excel software (Eq.3).As shown in this figure, the baseline scenario includes high degradation in some hexagonal units.However, the majority of its units experienced a degradation degree of less than 10 in the base year of 2015.For a more effective comparison, the classification of degradation degrees for both the baseline and future scenarios was categorized into four relative groups.This representation illustrates those scenarios 1, 2, and 3 are projected to experience increasing levels of degradation in the future.The degradation rate per hexagonal unit is labeled.
The comparison of the degradation rate in hexagonal units for each of the baseline and future scenarios is shown in Fig. 7.In general, the sum of degradation rate of hexagonal units from the base situation of 2015 with a value of 270.33 will reach 229.23, 250.60, and 1682.80 in scenarios 1, 2, and 3, respectively.

Patterns of landscape change
Examining the images of 1985-2000, 2000-2015, and 1985-2015 reveals significant changes in the land cover classes of Hoor Al-Azim wetland and its surroundings, especially in the water level.The decrease in the wetland's water level has resulted in changes in other covers, like dense and sparse vegetation.Throughout these periods, the wetland has undergone significant degradation, as evidenced by a marked reduction in water bodies from 207,511 ha in 1985 to 50,264 ha, indicating a decreasing trend.Additionally, the area of salt lands decreased from 33,181 ha in 1985 to 9628 ha in 2015 due to receding water and the growth of vegetation, such as reeds and cholan, along the wetland's shore.
Consistent with findings from studies such as those by Yu et al. 65 and Zhu & Gong 66 , other wetland regions have similarly experienced a reduction in water, with human activity identified as the main threat to these ecosystems.Conversely, the wetland area has increased in the Abitibi-Témiscamingue region of Canada, attributed to the minimal impact of human activities in this area 44 .
Various factors in the anthropogenic process have led to changes in the Hoor Al-Azim wetland's ecosystem.The main source of water for Hor comes from its upstream rivers, including the Karkheh in Iran and the Tigris and Euphrates in Iraq.Any interference, whether human or natural, with the flow rate of these rivers will significantly impact the water requirements of the wetland.In Iran, the Karkheh River and its branches supply most of the water to the Hoor Al-Azim wetland 39 .The construction of the Karkheh Dam in 1997 and its impoundment during 1998-2002 were of the most important factors in reducing the level of the wetland in such a way that the flow entering the wetland through the Karkheh River was 5450 million m 3 annually on average before the construction of the dam, which has decreased to 4700 million m 3 since 1998-2002.The evaporation rate of the wetland is 1.9 million m 3 per km 367 .Human factors, including resource threat and pollution since the war between Iran and Iraq, have degraded the Hoor Al-Azim wetland 68 .Factors like road construction, oil drilling, lack of limitations, and loss of effective protection have led to its drying and turn it into an area with dense, www.nature.com/scientificreports/sparse vegetation or salt land.The southern part has been affected by the Azadegan oil field development plan, while operations in the northern part are still ongoing 5 .Key strategies for wetland restoration include preventing habitat fragmentation by avoiding road and dike construction, harnessing local community support for wetland conservation and ecotourism, and utilizing national and international funding to provide protective infrastructure.Additionally, strategic zoning can optimize resource allocation for agriculture and ecotourism.

Dynamics modeling of wetland
The success and performance of simulation models based solely on quantity cannot provide researchers with information on spatial criteria and morphological features.Therefore, a spatially explicit scheme is needed to evaluate the performance of spatially explicit predictive models.In addition, this plan is expected to offer insight into the model's behavior.Therefore, this research proposed an advanced method based on the study of the landscape pattern to evaluate the model's performance 26 .According to the characteristics of the studied area, the hybrid approach used in this study was beneficial for comprehending the intricate dynamics of the landscape pattern of the Hoor Al-Azim wetland in the region.The MOLA algorithm used by LCM was successfully applied in the dynamic modeling of wetland changes.This study also evaluated the comparative simulation capabilities of integrated MLP_Markov and CA_Markov models in terms of ecology using the FOM index.Using this index is one of the most valuable and accurate methods in comparing and confirming the accuracy of simulation models, and it has been used in numerous studies.Based on the results of the accuracy evaluation, the FOM index value in the current study was higher for the MLP_Markov model than for the CA_Markov model.The rate of this index was found to be 14.73% in the comparison between the ground truth map and the simulated map of 2015 in Tianjin, China 69 .In the study of Zabihi et al. 70 , this index was reported to be 10.53% in the Talar watershed of Iran, which indicates that the accuracy of modeling in the present study is higher than those in the two mentioned studies.Kourosh Niya et al. 27 used an adding/deleting approach to increase the accuracy of hybrid models and assessed the accuracy of different models using the FOM index in Iran's Qeshm Island.The value of the FOM index when using this approach was 7.8, which was better than when not using it.Shoyama et al. 71 developed high-resolution land use maps by utilizing a GIS-based vegetation database.They also created scenarios to project potential land changes in Japan by 2050 on a national scale.The value of the FOM index in that study was 0.58%, which is significantly lower than its value in the present study.

Estimating the severity of degradation
The study generated different scenarios for 2030 and evaluated their impact on the landscape comparison to the baseline situation in 2015.By utilizing a degradation model and a network of hexagonal units, landscape metrics were used to estimate the different conditions resulting from each policy adopted.In a study conducted in East India, researchers found that landscape criteria and spatial characteristics of land use/land cover are effective tools for identifying land fragmentation and degradation 72 .Their analysis of the Chatra wetland revealed significant deterioration, with a 60% degradation over 28 years, due to urban construction.The transformation and change of natural land cover around the wetland have a destructive impact on the wetland ecosystems, mirroring the results of this research.Although, the education and social learning of people as stakeholders of land use change have been recognized as effective in improving resilience under drought conditions in Iran 73,74 .
While numerous studies have analyzed historical degradation rates in wetland ecosystems, there has been limited research projecting degradation into future scenarios.Hu et al. 8 , provided a framework for extracting the drivers of wetland degradation and found significant differences in degradation rates across scenarios, similar to the findings of this research.Evaluating wetland degradation has been a critical issue, and Cui et al. 6 proposed a framework that showed rapid ecological degradation in the Yancheng National Natural Reserve wetland in China over 30 years.The innovative approach used in this study, combining landscape metrics and land cover dynamics, offers a new method for evaluating future wetland degradation.The severity of degradation in management scenarios can be influenced by land managers, who can change metrics by altering land covers.Suitable land covers, such as water bodies, humid land, and dense vegetation (with increased values of CA and MPS metrics) contribute to lower degradation, while unsuitable land covers like salty lands and sparse vegetation (with decreased values of these metrics), lead to higher degradation.Also, less degradation caused by the decrease in ED, MPFD, and MSI values in suitable land covers and an increase in the value of these metrics in unsuitable land covers.On the other hand, the high percentage of water, humid land, dense vegetation and low percentage of salty lands and sparse vegetation in hexagonal units reduces the severity of degradation.
Awareness of different strategies and their consequences in various scenarios is crucial for understanding the future situation of a region.Territorial managers and custodians can avoid unfavorable paths and move towards a more favorable future by adjusting their strategies.Simulation results of scenario 1, an optimistic scenario, show that water conservation and preventing the conversion of water bodies improve or minimize degradation compared to the baseline situation.Preserving water requirements, preventing drilling and oil exploration, and dredging waterways are important considerations for achieving this goal.In the intermediate scenario, the rate of degradation also decreases.If the current historical scenario trend continues, the degradation situation will worsen significantly by 2030.Protective measures have been taken for sustainable wetland utilization and conservation in Iran and Iraq.Three units for cooperation and monitoring, education and awareness, and science and technology are working to implement the Ramsar Convention on Hoor Al-Azim Wetlands (Iran) and Horalhawizeh (Iraq).They tried to develop a mechanism for effective communication between the Center for Restoration of Marshes and Wetlands and oil companies to mitigate the impact on the wetland.
Change in Hoor Al-Azim wetland is a dynamic process, and uncertainties can affect its trajectory of wetland.The model was utilized under the assumption that the explanatory drivers of wetland loss will not vary significantly within the next 15 years 75 .The exclusion of highly dynamic variables such as roads may affect the outcomes of the model.It is proposed that new roads will be added to the model as dynamic input.For the utilization of precipitation and temperature data, maps at a 1 km resolution were employed, despite the lack of availability in the desired scale.This may introduce a degree of uncertainty in the model's input due to the scale of these maps.For future studies, it is necessary to examine the changes in government policies and economic development in the region that can affect the results of the transition potential modeling.The accuracy of the simulated land cover maps is strongly related to the accuracy of input data layers for the model, such as historical land cover maps.Temporal and spatial heterogeneity make wetlands hard to map with remote sensing imageries.The coarse classification of wetlands may also introduce bias, as Hoor Al-Azim wetland is a permanent freshwater wetland surrounded by reed that may overlap with sparse or dense vegetation classes.A finer and more accurate classification of wetlands is needed to capture them.The accuracy of classified land cover maps extracted from remotely sensed data with moderate spatial resolution of the Landsat images can lead to lower reliability than that of high-resolution satellite images in simulated maps.Therefore, it is suggested to use sensors with higher spatial resolution than Landsat images to produce more accurate land cover maps.The scientific findings of the current research can also be used as input in studies analyzing the impact of spatio-temporal future land cover scenarios on ecosystem health assessment.The study of determining the water requirement of the wetland under the impacts of different future land cover scenarios and providing optimal solutions for supplying the water demands of the wetland should be investigated in the future.

Conclusion
The present study aims to develop a highly accurate model for simulating future scenarios of land cover changes in Hoor Al-Azim International Wetland.The methodological approach employed involves a comparative analysis between the hybrid model-integrating an artificial neural network, Markov chain, and multi-objective optimization-and the cellular automata-Markov chain model.It was hypothesized that the LCM model would have a higher transition accuracy potential than the CA-MC model due to the use of artificial neural networks in modeling.The model-generated maps underwent successful validation via the FOM index application.The LCM model allowed us to generate three scenarios for 2030 concerning achieving higher accuracy.Scenario 1 (Water Conservation) will include the largest percentage of water bodies in the future.During 1985-2015, a decreasing trend was seen in the water bodies, humid lands, and dense vegetation, while sparse vegetation and salt lands demonstrated an increase in Hoor Al-Azim wetland.Evaluation of degradation in hexagonal units also showed the least degradation in the Water Conservation scenario compared with the other two scenarios in 2030.Therefore, appropriate and timely management measures and planning should be done by decision-makers to create sustainable development and protection of the wetland, along with environmental protection strategies, to reduce the intensity of changes.
We applied the Cramer's V coefficient, a chi-square-based metric of association between two variables, to quantify the significance of each variable in elucidating changes in the wetland.The results of Cramer's coefficient, which indicates the relationship between explanatory variables and land cover classes, showed that the DEM with a total Cramer's coefficient of 0.1923 had the least impact on the changes of Hoor Al-Azim wetland, while annual average temperature (0.3729) indicated the greatest impact on the changes of the region.Given that temperature and precipitation significantly influence the dynamics of the Hoor Al-Azim wetland, it is crucial to assess the effects of projected climate change scenarios and human-induced activities on the wetland's future state, with particular emphasis on water yield and water requirement.The accuracy rate of modeling results using artificial neural networks in different sub-models was also obtained between 0.64 and 0.96, demonstrating its high ability to model transition potential.Although the simulation results indicated a decreasing trend at the wetland water level, there are still more considerations for future modeling studies.For example, access to better spatial resolution in satellite imagery provides more detailed information about different class types that may not be detected by LANDSAT-based classification.
Rapid and extensive land cover changes in Hoor Al-Azim International Wetland might have serious environmental effects.For example, changes in the hydrologic regime and potential biological responses may cause changes in habitat quality and the abundance of native species and wintering birds in the wetland.To enhance the restoration and ongoing conservation of the wetland, practical measures could include curbing livestock overgrazing, raising awareness of the wetland's diverse benefits to ensure its protection, engaging local communities as environmental stewards, establishing expert-led monitoring and evaluation programs, and forming a protective buffer zone around the wetland by environmental agencies to thwart harmful human encroachment.Additionally, legal enforcement can be utilized to limit the spread of oil company operations and prevent the discharge of oil and urban wastewater into the wetland, coupled with the establishment of a crisis management center by Iran and Iraq.It is suggested to model water yield and sediment retention ecosystem services under the influence of land cover changes in subsequent research studies.Also, the change in the temperature and precipitation pattern under future greenhouse gas release scenarios should be modeled and added as a factor in the modeling of land cover changes.Also, using the results and method applied in this research can warn decision-makers and territorial managers about the process of anthropogenic changes and help them take action, plan, and sustainably protect the wetland to reduce the intensity of changes that lead to degradation.In general, the developed process can act as a decision-support system for managers and decision-makers.In this way, managers can understand the degradation effects of management plan alternatives in the units of study area before implementation with the help of the LD index calculation.Another use of our result in management is targeting a certain range of the degradation index according to the time and financial limitations of any management plan alternatives.

Figure 1 .
Figure 1.Geographical location of the study area.
Hit (the model predicts the change between 2000 and 2015 and the ground truth map also shows the change from 2000 to 2015), Miss (the model did not predict any change between 2000 and 2015 while the ground truth map has changed), and False Alarm (the model predicts the change during the years 2000 to 2015 while the ground truth map shows no change) are included.

Figure 4 .
Figure 4. Modeled land cover maps of 2015 using hybrid (a), CA-MC approaches (b), and Hit, Miss, and False Alarm areas in evaluating the accuracy of the modeled maps.

Table 1 .
Land cover classification scheme.

Table 2 .
63finition and characteristics of landscape metrics calculated in the research63.

Table 3 .
Coding of the intensity of degradation of land covers in the region according to the calculated landscape metrics.

Table 4 .
Results of Cramer's coefficients between variables and land cover classes.

Table 5 .
The matrix of resulting area transition obtained in the CA-MC approach.